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given. I focus on fluctuation measurements and on theoretical ideas which 
appeared experimentally fruitful. 
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1. Introduction 

With the advent of large acceptance detectors it became possible to ob- 
serve not one but tens or even hundreds of particles produced in a single col- 
lision of relativistic nuclei. Such a multi-particle state constitutes an event 
corresponding to a single high-energy collision. Event-by-event analysis is 
potentially a powerful technique to study relativistic heavy-ion collisions, as 
magnitude of fluctuations of various quantities around their mean values is 
controlled by system's dynamics. For example, the energy and multiplicity 
fluctuations of many body system are related to, respectively, the system's 
heat capacity and compressibility. The two susceptibilities strongly depend 
the system's state and they experience dramatic changes at phase transi- 
tions. So, measuring the fluctuations we can learn about effective degrees 
of freedom of the system and their interactions. 

Since the early 1990s the event-by-event physics has grown to a broad 
field of active research of relativistic heavy- ion collisions, see the review 
articles [U [2]. In the following I overview the achievements and failures; 
I discuss difficulties and future perspectives of the event-by-event physics. 
I mostly present experimental results and I focus on the theoretical ideas 
which appeared to be experimentally fruitful. Although I have tried to cover 
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the whole field, the choice of the results to be discussed is to some extend 
subjective. 



The first attractive idea of event-by-event analysis was formulated by 
Reinhardt Stock [3] who suggested to look for 'interesting' classes of events. 
The interesting events were meant the collisions where the quark-gluon 
plasma is produced or the collisions of exceptionally high multiplicity or 
energy density etc. Imagine there are 'hot' events with the temperature 
significantly higher than the average one. Let us further assume the event 
temperature can be quantified by M(pr) which is the transverse momentum 
averaged over particles from a given event. It is defined as 



where N is the event's multiplicity. If the 'hot' events indeed exist, then 
the distribution of M(j>t) should reveal it. Fig. 1 shows a typical example 
of the distribution of M{pt)- The measurement was performed in central 
Pb-Pb collisions at 158 AGeV by NA49 Collaboration at CERN SPS @|. 
As seen, no 'hot' events are observed - the M(pr) distribution is of boring 
Gaussian shape. 

Fig. 1 shows another typical feature of event-by-event distributions. 
Namely, the distribution of M(pt) obtained for the so-called 'mixed' events, 
where every particle is taken from a different event, is nearly identical with 
that obtained for real events. Since there are no inter-particle correlations 
in mixed events by construction, the similarity of the two distributions pre- 
sented in Fig. 1. shows that particles in real events are mostly independent 
from each other. The same is suggested by the Gaussian shape of the dis- 
tribution. The fluctuations present in mixed events are called statistical 
and the fluctuations, which remain after the statistical fluctuations are sub- 
tracted, are called dynamical. 

Since potentially interesting information encoded in dynamical fluctua- 
tions is not easily seen in the event-by-event distributions we have to use 
more subtle methods to infer it. So, in the two next sections I discuss 
quantities to be measured. 



2. Early Days Motivation 




»=i 



(1) 



3. Measurable Quantities 



In thermodynamics we have extensive quantities such as energy or parti- 
cle multiplicity, which are proportional to the system's volume, and intensive 
quantities such as temperature or various densities, which are independent 
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Fig. 1. The distribution of transverse 
momentum M(px) measured in cen- 
tral Pb-Pb collisions at 158 AGeV by 
NA49 Collaboration. The histogram 
and points correspond to, respectively, 
the mixed and real events. The figure 
is taken from [4] 
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Fig. 2. The multiplicity distribution 
measured in Pb-Pb collisions at 158 
AGeV by WA98 Collaboration. The 
minimum bias data and three classes of 
central events are shown. The figure is 
taken from [5] 



to the system's size. One is tempted to introduce analogous quantities in 
event-by-event physics of relativistic heavy-ion collisions. The number of 
participants is a natural measure of the size of system which emerges in 
heavy-ion collisions. Then, the quantities like the energy carried by all pro- 
duced particles or particle multiplicity are approximately proportional to 
the number of participants and thus they are extensive. 

Fig. 2 shows the distribution of multiplicity of charged particles produced 
in Pb-Pb collisions at 158 AGeV at different centralities. The measurement 
was performed by WA98 Collaboration at CERN SPS [5]. The collision 
centrality is defined as a percentage of total inelastic cross section <j mel of 
nucleus- nucleus collision. The centrality of n% corresponds, roughly speak- 
ing, to the collisions with impact parameters from such an interval [0, b] that 
7t6 2 is n% of a 1 . As seen in Fig. 2, the smaller centrality (more central 
collisions) , the higher average multiplicity and the smaller width of the dis- 
tribution. The most upper curve corresponds to the minimum bias events 
when the collisions are collected with no selection - there is no centrality 
trigger condition. Fig. 2 shows that measurements of extensive quantity like 
multiplicity are not very informative, as the results crucially depend on a 
trigger condition. 

The quantities like M(px) are expected to be analogous to thermody- 
namic intensive quantities. Fig. 3 shows the distribution of transverse mo- 
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Fig. 3. The distribution of transverse momentum M(j>t) measured in four central- 
ity classes of Au-Au collisions at y / J/\w = 130 GeV by PHENIX Collaboration at 
RHIC. The figure is taken from [5J. 



mentum M(px) measured in central Au-Au collisions at ^Jsnn = 130 GeV 
by PHENIX Collaboration at RHIC [6j. As seen in Fig. 3, the average 
value of M(pt) is indeed approximately independent of the system's size 
(centrality) but the width of the M(px) distribution clearly depends on the 
system's size. And it is unclear whether the width simply depends on the 
trigger condition or it results from dynamics of nuclear collisions. 



4. Fluctuation measures 

In light of previous considerations it is desirable to construct a fluctu- 
ation measure which is truly intensive and it vanishes in absence of inter- 
particle correlations. Several quantities, which satisfy these conditions, have 
been proposed but I focus on the measure <I> introduced in [7]. It is con- 
structed as follows. One defines the single-particle variable z = x — x with 
the overline denoting averaging over a single particle inclusive distribution 
which is performed as 
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where iVj. is the particle multiplicity in k— th event, M is the number of 
events and A^totai is the total number of particles in M events. Thus, we 
sum over events and over particles from every event. The event variable Z, 
which is a multiparticle analog of z, is defined as 

N 

z = J2(xi -x), (3) 
i=l 

where the sum runs over particles from a given event. The averaging over 
events is 

(^> = ^E Z *- ( 4 ) 
fc=i 

One observes that by construction (Z) = 0. Finally, the measure <I> is 
defined in the following way 

The measure $ possesses two important properties: 

• when particles are independent from each other - there are no corre- 
lations among particles coming from the same event, the $— measure 
vanishes identically; 

• when particles are emitted by a number of identical sources, which 
are independent from each other, $ has the same value as for a single 
source independently of the distribution of the number of sources (<£ 
is strictly intensive). 

Due to the first property $ is exactly zero for mixed events. Because 
of the second property it is strictly independent of centrality in a broad 
class of models of nucleus-nucleus collisions where produced particles origi- 
nate form independent sources. The models include the Wounded Nucleon 
Model [8] and various models where a nucleus-nucleus collision is treated as 
a superposition of independent nucleon-nucleon interactions. In more real- 
istic transport models like HIJING [9], VENUS [10], UrQMD [llj or HSD 
|12j , there is an admixture of secondary interactions which break down inde- 
pendence of nucleon-nucleon interactions. However, $ is still approximately 
independent of centrality within these models. 

As already mentioned, several other fluctuation measures were intro- 
duced. In Ref. [13], see also [H] , it was proposed to use 

a% Q = ((X - (X)) 2 ) - j^T^W , (6) 
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where X is the event variable 

i=l 

The authors of [TS] advocated the measure 



1 N 



a 



2 

dyn 



S = sgn(^ y J^^. (8) 

We also mention here the quantity F introduced in [16] which is defined in 
the following way. One obtains the scaled dispersion 



for real events and for mixed events, and then one computes 

p ^data k-'mixed (10) 

k-'mixed 

The fluctuation measures cij yn , S and F similarly to $ vanish in the 
absence of inter-particle correlations. However, none of these measures is 
strictly intensive as $ is. Knowing the average multiplicity (N), the mea- 
sures <J>, c| , £ and F can be approximately recalculated one into another. 



5. Transverse Momentum Fluctuations at SPS 

Transverse momentum fluctuations in nucleus-nucleus collisions at SPS 
energies were measured by NA49 |17j and CERES [15] Collaborations. 
Fig. 4, which is taken from [17], shows the data on p-p, C-C, Si-Si and 
Pb-Pb collisions at 158 AGeV. The fluctuations are measured by means of 
at various centralities determined by number of wounded nucleons. In 
Fig. 1 we can hardly see the difference between the real and mixed events, 
Fig. 4 clearly demonstrates presence of dynamical fluctuations and thus 
it proves sensitivity of the measure. The magnitude of the dynamical 
correlations is quite small (<I> < 8 MeV) when compared to the dispersion 
of the inclusive transverse momentum distribution (the second term in the 
definition (0)) which varies within 200 MeV < a PT < 250 MeV [T7]. So, the 
dynamical correlation is a few percent effect. 

We observe that the fluctuations are different for positive and negative 
particles. It is not surprising as the negative particles are nearly all negative 
pions while the positive particles include sizeable fraction of protons (the 
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Fig. 4. Q(pt) as a function of wounded 
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Fig. 5. Two-particle correlation plot of 
the cumulant variables x%, X2 in cen- 
tral Pb-Pb collisions at 158 AGeV. The 
figure is taken from [15] . 



measurement shown in Fig. 4 was performed in the forward hemisphere). 
We also observe the centrality dependence of $ with the maximum at rather 
peripheral collisions. 

Although, the measure $ is sensitive to various dynamical fluctuations, 
one needs more differential observables to identify a nature of the fluctu- 
ations. For such a purpose one can use the two-dimensional plot of the 
cumulant variables x±, X2 proposed in [13]. Following [18], one defines the 
cumulat variable 

rpr 

x(pr) = / dp' T P(p' T ). (11) 
Jo 

where P{pr) is the inclusive distribution of pr- Since P{pr) is normalized 
to unity, < x < 1. And now one finds a point {x\,X2) for every pair of 
particles from the same event and constructs a two-dimensional plot such as 
shown in Fig. 5 [19] . In the absence of any correlations the plot is flat and 
various correlations generate different patterns in the plot. The example 
shown in Fig. 5 proves an existence of positive correlation among particles 
of the same px which is signaled by the ridge along the diagonal. Obviously 
the correlation is due to the Bose-Einstein statistics of identical pions. 

The results of even more differential analysis performed by CERES Col- 
laboration [20] are shown in Fig. 6. The pairs of particles, which contribute 
to the correlation plot, are divided into classes according to the relative 
azimuthal separation of the two particles A<I>. As seen in Fig. 6, the pat- 
tern of correlation qualitatively changes with A<£. For the small separation 
0° < A$ < 30° we observe the Bose-Einstein correlation, but for the max- 
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Fig. 6. Two-particle correlation plot of the cumulant variables x\, X2 in central Pb- 
Au collisions at 158 AGeV. The left and right figures, which are taken from [2D], 
correspond to the relative azimuthal separation of the two particles 0° < A$ < 30° 
and 150° < A$ < 180°, respectively. 

imal separation 150° < A<I> < 180° the correlation is presumably caused 
by the event-by-event fluctuations of the slope of transverse momentum 
distribution. 

The correlation plots shown in Figs. 5, 6 are indeed informative but 
still there is a correlation which is not clearly seen in these plots. This is 
the correlation of the event's transverse momentum and event's multiplicity 
which was observed long ago in p-p collisions at 205 GeV |21| . The correla- 
tion appears to be sufficiently strong to give a significant, if not dominant, 
contribution to shown in Fig. 4 |22j . 

We conclude this section by saying that the dynamical transverse mo- 
mentum fluctuations in heavy-ion collisions at SPS are of various physical 
origin but their total magnitude is quite small. 

6. Transverse Momentum Fluctuations at RHIC 

Transverse momentum fluctuations in nucleus-nucleus collisions at RHIC 
were measured by PHENIX Collaboration [T5] using F, see Eq. (|1Q|) . and by 
STAR Collaboration [23] using <?\ jn , see Eq. ([6]). Fig. 7, which is taken from 
[16] . shows the centrality dependence of pr fluctuations which appears to 
be similar to that at SPS. The magnitude of the fluctuations is bigger. The 
measurement performed by STAR Collaboration [23J, which can be easily 
recalculated into $(pt)> shows that 3>(pr) exceeds 50 or even 70 MeV at 
top RHIC energy. However, it is difficult to quantitatively compare results 
from different experiments because the measured fluctuations depend on the 
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Fig. 7. F(px) as a function of number 
of participating nucleons in Au-Au col- 
lisions at y / s]v7V = 200 GeV. The figure 
is taken from [TBI. 
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Fig. 8. F(px) as a function of upper pt 
cut-off for Afpart » 150 in Au-Au colli- 
sions at ^/snn — 200 GeV. The figure 
is taken from [TBI. 



acceptance which differs from experiment to experiment. 

It was observed in [16] that the px fluctuations are dominated by parti- 
cles with relatively high pr- Fig. 8 shows F(pr) as a function of upper pt 
cut-off for the centrality corresponding to the maximal fluctuations. For a 
given p™ ax only particles with p? < p™ ax are taken into account. As seen, 
F(pt) grows fast with p™ ax and consequently it was claimed [16] that the px 
fluctuations are due to jets or mini-jets. The claim, however, was questioned 
in [21] where it was argued that the data from Fig. 8 can be reproduced 
within a statistical model with multiple clusters or fireballs which move at 
some collective velocities, correlating the momenta of particles belonging to 
the same cluster. Thus, similarly to the situation at SPS, there is no unique 
interpretation of dynamical pr fluctuations at RHIC. 



7. Thermodynamic fluctuations 

As mentioned in the Introduction, fluctuations in many body systems 
carry information about the system's state and its dynamics. Assuming that 
the strongly interacting matter produced in relativistic heavy-ion collisions 
is in thermodynamic equilibrium, it was suggested |26^ [27] to measure the 
temperature fluctuations. Then, using the relation 



(12) 
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which is discussed by Landau and Lifshitz [25], one can infer the system's 
heat capacity at fixed volume V and particle number N 

c^m , (.3) 



V.N 



where U is the system's energy. The relation (|12|) . however, is actually 
very controversial [28[ [29] and its status is rather unclear. Not entering 
the details, I think that the relation (|12p cannot be used, as long as the 
thermometer to measure the temperature fluctuations is not specified [30] . 
A similar idea [31] was to infer the compressibility 



dp 
dV 



(14) 

T,(N) 



where p is the pressure, from the multiplicity fluctuations due to the relation 

m 

(N 2 ) - (AO 2 = (15) 

An experimental problem here is to measure the multiplicity fluctuations at 
fixed system's volume. 

Only the third idea to study electric charge fluctuations in relativistic 
heavy-ion collisions appeared to be experimentally relevant. The fluctua- 
tions are related to the electric charge susceptibility [2] as 

(Q 2 ) - (Q) 2 = TV XQ , (16) 

with 

« s -(^L- (17) 

where F is the free energy and \iq is the chemical potential responsible for 
the electric charge conservation. Eqs. (|16l [TTl) do not look very exciting at 
first glance but it was sharply observed [32j [33] that the susceptibility (fT7|) 
is very different in the quark phase and in the hadron one. 

To explain this statement, let me consider the classical ideal gas of par- 
ticles of chargers ±q (measured in the units of elementary charge). The 
system's charge is then Q = q(N + — A_). We introduce 5Q = Q — (Q) and 
5N± = N± — (N±) and we compute the charge fluctuations as 

((5Q) 2 ) = Q 2 m + - 5N_) 2 ) = q 2 (((5N + ) 2 ) + (5N_) 2 ) - 2(5N + 5N_)). 

Since in the ideal classical gas ((SN±) 2 } = (N±) and (5N+5N-) = 0, one 
finds 

{{5Q)2) - a 2 (18) 

1jv> — 9 ' (18) 
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Fig. 9. Electric charge fluctuations 
quantified by <& 9 as a function of rel- 
ative charge multiplicity in central Pb- 
Pb collisions at SPS for several collision 
energies. The figure is taken from [34] 
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Fig. 10. Electric charge fluctuations 
quantified by A$ g = <i> g - $^ cc as a 
function of relative charge multiplicity 
in central Pb-Pb collisions at SPS for 
several collision energies. The figure is 
taken from 134 . 



where (N) = (N + ) + (iV_). As seen in Eq. (fl~8j) . the charge fluctuation per 
particle equals the particle's charge squared. 

One easily derives the formula analogous to Eq. (|18p for the ideal clas- 
sical gas of pions composed of 7r + , tt~ , tt° and for the quark-gluon plasma 
being a mixture of ideal classical gases of quarks of different charges and of 
neutral gluons. Using the system's entropy S instead of the total particle 
multiplicity (N), one finds [2] 

((SQ) 2 ) _ { I for pions, 
S ~ \ i for QGP. 

It was argued in [321 [33j that the charge fluctuations generated in the quark 
phase are frozen due the system's fast hydrodynamic expansion and that the 
entropy, which is mostly produced at the very early, preequilibrium stage 
of the collision, is approximately conserved during the hydrodynamic evolu- 
tion. Then, a measurement of the ratio (|19p should clearly show whether the 
quark-gluon plasma is produced at the early stage of relativistic heavy-ion 
collisions. 

Soon later the electric charge fluctuation were measured experimentally. 
Figs. 9 and 10 show the results obtained at SPS by NA49 Collaboration [M] 
using the measure $ defined by Eq. ([5]). (-/V c h)tot and (-/V c h) are the average 
charge particle multiplicities in, respectively, the full (47r) acceptance and 
in a given phase-space domain under study. As seen, the results, which are 
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Fig. 11. Electric charge fluctuations quantified by i/J 5 ^ 1 as a function of pseudora- 
pidity density of charged particles in nucleus-nucleus collisions at RHIC collision 
energies. The vertical axis shows multiplied by the pseudorapidity density of 
charged particles. The figure is taken from |39j . 



essentially independent of the collision energy, follow the trend dictated by 
the global charge conservation (GCC) corresponding to 



(-Nch)toi 



$jcc = / : _ _A^^ _ 1 (20) 



The formula (f2Uj) derived in [35] is actually approximate as it is derived 
under the assumption that the total system's charge Z vanishes. It is, 
however, a reasonable approximation of the exact formula derived in [36] 
when Z < (Ach)tot- 

Fig. 10 shows the electric charge fluctuations when the effect of the 
global charge conservation is subtracted that is there is presented A$ g = 
& q — $^ cc . Fig. 10 also shows the levels of charge fluctuations in the quark- 
gluon plasma and in the hadronized system, both computed in a rather 
simplified model. As seen, the observed fluctuations agree very well with 
the hadron gas prediction. 

The electric charge fluctuations were measured at RHIC by PHENIX 
[37] and STAR [381 EH] Collaborations. The results shown in Fig. 11, which 
is taken from [39J, are rather similar to those obtained at SPS. However, the 
STAR Collaboration used the measure z/J 5 ^ 1 to quantify the electric charge 
fluctuations. It is defined as 

dyn _ (N + (N + - 1)) (iV_(jV_ - 1)) (N + N_) 

+-- (iV + )2 + (AT_)2 <JV + ><J\r_>- 1 ] 
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is sensitive only to the dynamic fluctuations in this sense that it van- 
ishes when the fluctuations of both N + and iV_ are Poissonian. 

As seen in Fig. 11, the observed fluctuations are not only bigger than 
those in QGP but they are even bigger than those in the hardon resonance 
gas. Although we have good reason to claim that the quark-gluon plasma is 
produced at the early stage of relativistic heavy-ion collisions at RHIC, the 
final state charge fluctuations do not signal the presence of the QGP. Most 
probably the fluctuations generated at the plasma phase are simply washed 
out during the subsequent system's evolution. The fact that the observed 
charge fluctuations are bigger than those in the hardon resonance gas is pre- 
sumably caused by a relatively small acceptance of the measurement. When 
a significant fraction of particles originate from neutral resonances, which 
decay into one positive and one negative particles, the charge fluctuations 
are reduced, when compared to the Poissonian fluctuations, if both particles 
from the decay are observed [35J. When the experimental acceptance is so 
small that typically only one particle from a resonance decay is registered, 
the the electric charge fluctuations remain Poissonian. 



8. Balance functions 

In the previous section I discussed bulk fluctuations of electric charge 
which at the end appeared to be not very informative. Here I am going to 
present a very interesting idea [4Q|, [4"T] to measure correlations of the electric 
charges in rapidity by means of the so-called balance functions defined as 

(JV + _(Ay))-(iy_-(Aj/)) (N_ + (Ay)} - (N ++ (Ay)) l 

(JV_(Ay)) + (N+(Ay)) \ { > 

where (N±(Ay)) and (N±±(Ay)) are, respectively, the average number of 
positive or negative particles and the average number of pairs of particles 
of given charges within the rapidity (or pseudorapidity) interval Ay (Ary). 
The balance functions were argued [401 [42] be sensitive to a hadronization 
mechanism. The width of the balance functions was expected to be bigger, 
when the hadronization proceeds via the break-up of strings as in p-p colli- 
sions, than when the quark-gluon plasma hadronizes due to the coalescence 
of constituent quarks. 

The balance functions were measured in Au+Au collisions at RHIC |43j 
and in Pb-Pb collisions at SPS [H], see Fig. 12 and 13. The balance func- 
tions for peripheral collisions appeared to have widths consistent with model 
predictions based on a superposition of nucleon-nucleon scattering. Widths 
in central collisions were smaller, consistent with trends predicted by mod- 
els incorporating late hadronization due to the coalescence mechanism. Un- 
fortunately, the interpretation appeared to be not unique as the balance 
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Fig. 12. The balance functions in cen- 
tral and peripheral Au-Au collisions 
^/snn — 130 GeV. The figure is taken 
from [43]. 
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Fig. 13. The balance functions in Pb-Pb 
collisions at different centralities at 158 
AGeV. The centrality class 'Veto V cor- 
responds to the most central collisions. 
The figure is taken from [44] , 



functions were shown to be influenced by various factors [45 } 146} 147] . In 
particular, it was observed that the variation of the amount of transverse 
flow with collision centrality can reproduce [Uj the experimentally observed 
narrowing of the balance functions for central collisions. 



9. Multiplicity fluctuations 

As discussed in Sec. [3l the multiplicity measurements like that one pre- 
sented in Fig. 2 are not very useful, as the results crucially depend on the 
collision centrality. The situation is changed if the centrality condition does 
not result from specific features of a detector used in the measurement but 
if the centrality condition corresponds to a well denned physical criterion. 
Such measurements were performed by the NA49 Collaboration [IS] with the 
help of zero degree calorimeter which allowed one to determine the number 
of participating nucleons from a projectile (Np™l ectllc ) in a given nucleus- 
nucleus collisions. Fig. 14 shows the scaled variance (((N — (N)) 2 }/(N)) 
as a function of N^l°i ectlle in p-p and Pb-Pb collisions at 158 AGeV [4"5] . 
We observe a non-monotonic behavior of ((N — (N)) 2 ) / (N) which contra- 
dicts commonly applied models. In the Wounded Nucleon Model [8J, where 
produced particles come from wounded nucleons, which are assumed to be 
independent from each other, the scale variance is exactly independent of 

projectile _ Ag seen [q pig _ ^ ^ transport mode l s HIJING 0, VENUS 

[TO] . UrQMD [TTJ or HSD [12] predict the approximate independence. It 
should be noted here that although the scaled variance is a non-monotonic 
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Fig. 14. The scaled variance of multiplicity distribution of negative (upper panel), 
positive (middle panel) and charged (lower panel) particles as a function of number 
of projectile participants in nucleus-nucleus collisions at 158 AGeV. The predictions 
of HIJING, VENUS, UrQMD and HSD models are also shown. The figure is taken 

from gg. 



function of A^° t j ectlle , the average multiplicity is simply proportional to 

Afpa°t° Ctlle [48j in agreement with the models mentioned above. Although 
there were several theoretical attempts [50J EH [52] to explain the data 
shown in Fig. 14, in my opinion, there is no reliable explanation. 

The multiplicity distribution at the most central collisions reveals an 
interesting feature. As shown in Fig. 15 taken form [53] it is narrower not 
only than the Poisson distribution but it is narrower than the multiplicity 
distribution obtained in the statistical model [54J which uses the Canonical 
Ensemble where the electric charge is exactly conserved. As seen in Fig. 16, 
this feature persists in a broad range of collision energies. Fig. 16 also shows 
that within statistical models one has to refer to a microcanonical ensemble 
to reproduce the scaled variance of multiplicity distribution. 

The multiplicity distributions discussed here appear to be associated 
with the transverse momentum fluctuations discussed in Sec. [6j As seen in 
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Fig. 15. The multiplicity distribution 
of negative charge particles produced 
in the most central Pb-Pb collisions at 
158 AGeV. The distribution is divided 
by the Poisson distribution of the same 
mean. The predictions of statistical 
models based on the Grand Canoni- 
cal and Canonical Ensembles are also 
shown. The figure is taken from [53j . 



0.95 



0.90 



0.85 




Fig. 16. The scaled variance of mul- 
tiplicity distribution of negative parti- 
cles produced in the most central Pb- 
Pb collisions as a function of collision 
energy. The predictions of statistical 
models based on the Grand Canonical, 
Canonical and Microcanonical Ensem- 
bles are also shown. The figure is taken 
from 1531. 



the definition of the measure <E> (|5"]), it depends on the multiplicity distri- 
bution. It was shown in [22] that the correlation of the event's transverse 
momentum and multiplicity, which is observed in p-p collisions [2TJ, com- 
bined with the non-monotonic scaled variance of multiplicity distribution 
shown in Fig. 14 approximately reproduces the pr fluctuations shown Fig. 4. 
Therefore, the similarity of Figs. 4, 14 is far not superficial. 



10. Elliptic Flow Fluctuations 

The elliptic flow is caused by an azimuthally asymmetric shape of the 
initial interaction zone of colliding nuclei. Consequently, it is mostly gener- 
ated in the collision early stage. Fluctuations of the elliptic flow were argued 
to carry information on very early stages of relativistic heavy-ion collisions 
[55] [56] . Large fluctuations of the elliptic flow were indeed observed at 
RHIC by PHOBOS [57J and STAR [58] Collaborations. However, STAR 
Collaboration claimed later on [59] that the magnitude of the fluctuations 
should be taken only as an upper limit due to the difficulties to disentangle 
the elliptic flow fluctuations and the contributions which are not correlated 
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Fig. 17. The relative fluctuations of 
the elliptic flow in Au-Au collisions at 
y/s NN = 200 GeV. The figure is taken 
from [57]. 




X 
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Fig. 18. Positions of wounded nucleons 
in the plane transverse to the beam in 
the Au-Au collision. The figure is taken 
from [62 . 



with the reaction plane. PHOBOS Collaboration has not retracted the data 
|57| . The whole problem is discussed in detail in the very recent review [BD] . 

As seen in Fig. 17, the relative v% fluctuations measured by PHOBOS 
Collaboration [57] are as large as about 40%. It appears, however, that the 
effect is dominated not by the dynamics but by simple geometrical fluctua- 
tions of the eccentricity of the interaction zone as suggested in [61]. Since 
the positions of nucleon-nucleon interactions fluctuate within the overlap 
region of the colliding nuclei as illustrated in Fig. 18 taken from [62], the 
eccentricity of the region fluctuates as well. Since the elliptic flow is pro- 
portional to the eccentricity, the relative eccentricity fluctuations directly 
contribute to the relative elliptic flow fluctuations. The calculations of the 
eccentricity fluctuations reproduce well the experimentally observed elliptic 
flow fluctuations, see e.g. [62]. Therefore, the hydrodynamic evolution of 
the system, when the elliptic flow is generated, seems to be fully determin- 
istic. The result is rather paradoxical if one remembers that the elliptic flow 
is mostly generated at a very early stage of the collision when the produced 
matter is presumably not in a complete equilibrium yet. 



11. Conclusions and Outlook 

A big volume of experimental data on event-by-event fluctuations in 
relativistic heavy-ion collisions has been collected for last fifteen years. Some 
results are indeed very interesting but the observed fluctuations are usually 
dominated by statistical noise as convincingly illustrated by similarity of 
mixed and real events. Theoretical expectations of large fluctuations cased 
by, say, phase transitions appeared to be far too optimistic but measuring 
of the fluctuations has also appeared rather difficult. 
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When single particle distributions are measured a detector inefficiency 
is not a serious obstacle. A number of undetected particles should be esti- 
mated and the single particle distribution is then easily corrected. In the 
case of correlation measurements, the effect of lost particles on the mea- 
sured correlation depends on how the lost particles are correlated with the 
detected one. Since the correlation is a priori not known, it is unclear how 
the observed correlation should be corrected. For this reason, the correlation 
measurements were usually performed in rather small acceptances where the 
detector efficiency is almost perfect. Then, the observed correlation signal 
does not need a correction for lost particles. However, dynamical correla- 
tions are usually strongly diluted due to a small acceptance. As an example, 
let me consider a multiplicity distribution. If we detect only a small fraction 
p of all particles, the observed multiplicity tends to the Poisson distribution 
when p — > 0. Consequently, we observe the Poisson distribution in a small 
acceptance independently of the actual distribution. I note that currently 
no more than 20% but typically only a few percent of all produced particles 
are used in event-by-event studies. 

Another problem of the current experiments is that the actual collid- 
ing system is not well known as an averaging over a centrality interval is 
performed. Such an averaging dilutes a potential signal, as most of charac- 
teristics of heavy-ion collisions strongly depends on centrality. Sometimes 
the centrality is estimated using produced particles which are analyzed. 
Then, the effect of autocorrelation has to be additionally removed from the 
data. 

The analysis of multiplicity clearly shows how important is a good de- 
termination of centrality. The multiplicity measurement presented in Fig. 2 
badly depends on experimental condition and thus is not very useful. When 
the collision centrality is so precisely measured that the number of partic- 
ipating nucleons from a projectile is known, the multiplicity distribution 
appeared to conceal very interesting features displayed in Fig. 14, 15. 

As the observed dynamical fluctuations are usually small, it is difficult to 
extract physically interesting information, it is even more difficult to work- 
out a unique interpretation. New theoretical ideas and reliable models are 
certainly needed but what the event-by-event physics really requires is, in 
my opinion, a new generation of experiments which will fulfill two impor- 
tant conditions: i) the acceptance is a sizeable fraction of 4tt, ii) the collision 
centrality is measured up to single nucleons participating in a collision. The 
future NA61/SHINE program at SPS is hoped to satisfy the requirements 
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